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Abstract 

Tentative modelling of the exploitation System of the ponds of Dombes (Ain, France). 

— The Dombes is a région in the Ain department, in the centre of which are many arti- 
ficial large ponds. These ponds are alternately exploited dry or flooded. The water 
phase (’évolage’) lasts from two to six years and is mostly devoted to the raising of 
carp, tench, pike and whitefish. The yearly fishing is accomplished by emptying the 
pond. The dry phase (’assec’) lasts one or two years and is devoted to the culture of oats, 
wheat, or maize. Two types of models are presented. The first one is an eco-energetical 
model, which describes the energy fluxes of the System in a stationary State. The energeti- 
cal cost of the food production of this ecosystem can be estimated. It is shown that the 
energetical efïiciency in this System is higher than in normal agricultural Systems. The 
second model is of a dynamic type. It describes the temporal variations of the System 
variables during the water phase. Six compartments are considered: phytoplankton, 
zooplankton, fish, détritus, sédiment, and minerai nitrogen. The dynamics of these six 
State variables is governed by a System of differential équations. The analytical form 
of the interactions, and the values of the parameters are taken from the literature. 
After simulating seven years, it is posible to compare the prédictions of the model 
to the real productions. The match is not perfect but qualitatively good. The model 
correctly predicts a graduai decrease of the yield during the water phase. 
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INTRODUCTION 

La Dombes est un plateau d’origine glaciaire situé dans le département de l’Ain 
(France). Au centre de ce plateau se trouve une région de 900 km 2 occupée par de nom¬ 
breux étangs artificiels. Ces étangs, exploités alternativement en eau et en assec, pro¬ 
duisent selon les années du poisson ou des céréales. L’objet de ce travail est la modéli- } 
sation de ce système d’exploitation agro-piscicole, plus particulièrement sous l’angle 
de la production alimentaire, afin d’en décrire le fonctionnement et d’en simuler la dyna¬ 
mique. 

D’une altitude moyenne de 280 m, cette région possède un climat de type rhoda- ( 
nien à nuance continentale. A Marlieux, la température moyenne pour la période 1964- 
1977 est de 19° C en juillet (maximum absolu 36,9° C) et de 2° C en janvier (minimum 
absolu — 23° C). L’insolation annuelle moyenne pour cette période est de 2000 h. Le 
régime des précipitations pour la période 1972-1982 atteint la moyenne annuelle de 
883 mm. 

La physionomie des étangs dombistes varie beaucoup selon leur végétation, leur 
profondeur et leur mode d’exploitation. Leur superficie oscille entre quelques hectares 
et une centaine. La profondeur est de 1 à 2 m vers la tête, à l’emplacement de la chaussée 
(digue), et de quelques centimètres à la queue et vers les rives. 

Aucun cours d’eau important ne traversant la région, les étangs sont alimentés 
par un réseau de fossés qui canalisent les eaux de ruissellement provenant des terres 
cultivées du bassin versant et les eaux de vidange des étangs supérieurs et mitoyens. 
Lorsque la pente est importante, les étangs sont disposés en chapelet dans l’alignement 
d’un thalweg. Une seule chaussée sépare la tête de l’étang supérieur de la queue de 
l’étang inférieur qui doit alors accepter les eaux sans pouvoir les détourner. 

L’origine des étangs est très ancienne puisqu’on en trouve la trace dans les textes 
dès le 13 e siècle. Rédigée en 1524, la « Coutume de Villars » précise un certain nombre 
de règles ayant trait à la création et au fonctionnement des étangs. Ce texte sera repris 
et complété par différents auteurs et en particulier Rivoire & Truchelut (1881) qui 
expliquent ainsi les conditions dans lesquelles ce système a pu être mis en place: «Il 
était loisible à tout le monde de construire un étang sur son terrain. Mais le propriétaire 
de la chaussée devait traiter avec ses voisins dont il inondait le terrain, soit en achetant 
leur fonds, soit en faisant des échanges, mais plus généralement, en payant les deux 
tiers de la valeur du terrain, et se soumettant d’abandonner ce même sol, la troisième 
année, aux propriétaires, pour le semer à leur profit. » 


FONCTIONNEMENT DU SYSTÈME D’EXPLOITATION DES ÉTANGS 

| 

Chaque étang possède un fossé de vidange et généralement un ou plusieurs fossés 
d’arrivée d’eau. Dans de nombreux cas celle-ci peut être détournée par une rivière de 
ceinture. Les eaux superflues sont évacuées par un trop-plein: l’ébie. Aucun étang ne 
pourrait se remplir ni se vider sans l’existence de ces multiples fossés qui communiquent 
entre eux, récupèrent les eaux de ruissellement et forment, dans chaque commune, un 
réseau hydrographique d’une complexité variable selon sa topographie, l’importance 
de la pente, le nombre et la taille des étangs. 

Lors de l’évolage, un nombre important d’étangs reçoit du fumier ou du lisier en 
automne (500 à 1000 kg/ha). Certains exploitants utilisent des engrais phosphatés 
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Fïg. 1. — Cycles et options possibles d’assolement pratiqués dans les étangs, exploités successivement en eau et en assec. L’évolage 
(période d’exploitation en eau) s’étend en moyenne sur quatre ans, avec une pêche chaque annee. L assec lui succédé pour une durée 

d’un à deux ans. 
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(300 kg/ha) ou des superphosphates (150 kg/ha). La teneur du sol en calcaire est voisine 
de zéro. En juin, les étangs peuvent recevoir du calcaire broyé (1 à 2 t/ha). 

Les étangs produisent des carpes ( Cyprinus carpio L.) et des tanches ( Tinca tinca L.) 
qui sont commercialisées après 3 ans, ainsi que du blanc (rotengles Scardinius erythro- 
phtalmus L. et gardons Leuciscus rutilus L.) et des brochets (Esox lucius L.). Les exploi¬ 
tants élèvent eux-mêmes les alevins de carpes dans de petits étangs peu profonds. 

Les chiffres d’empoissonnage sont établis avec précision et dépendent des caractéris¬ 
tiques du fond. Les quantités introduites à l’hectare sont en moyenne de 100 à 150 panots 
(carpes de 2 étés), 10 à 20 kg de tanchons, 15 à 30 kg de blanc et 10 à 20 têtes de brochets. 
Les alevins de tanches, blanc et brochets sont prélevés en cours de pêche et immédiate¬ 
ment reversés dans l’étang d’exploitation. Les panots, provenant d’étangs d’alevinage, 
sont introduits plus tard dans l’année. 

La pêche annuelle se fait par vidange complète de l’étang. Le moment venu, l’exploi¬ 
tant ouvre le thou (vanne) et vide partiellement l’étang de manière à ce que seul le bief 
(fossé creusé dans la partie basse de l’étang) soit encore en eau. Les poissons s’y re¬ 
groupent. Les pêcheurs groupés en deux équipes de part et d’autre du bief tirent les 
filets de manière à rabattre le poisson vers le bord. En moins de quatre heures, une 
équipe d’une vingtaine d’hommes peut pêcher, trier et peser les tonnes de poissons que 
fournit l’étang. Les rendements habituels des étangs dombistes sont de l’ordre de 200 
à 250 kg par hectare. 

Après la période d’évolage de quatre à six ans en moyenne, l’étang est mis à sec 
pour être cultivé. Il est alors labouré avant que le fond ne forme une couche trop dure. 
Lors de l’assec, le bief est dégagé car il s’envase rapidement durant l’évolage. Les 
cultures les plus fréquentes sont le maïs et l’avoine. La culture du blé, plus rare, ne peut 
se faire qu’en deuxième année d’assec, cette céréale devant être semée en automne. 
Les rendements en avoine peuvent atteindre 25 quintaux à l’hectare au cours d’une 
bonne année et être pratiquement nuis si les conditions sont mauvaises. En maïs, ils 
varient entre 50 et 80 quintaux à l’hectare, et en blé entre 25 et 50. 

La figure 1 schématise le fonctionnement du système intégrant les phases d’évolage 
et d’assec, avec les différentes possibilités d’assolement pratiquées dans les étangs. 

Dans ce travail, deux types de modèles sont présentés. Le premier est un modèle 
éco-énergétique qui décrit le fonctionnement de l’ensemble du système dans un état 
stationnaire sous l’angle particulier des flux d’énergie. Le second est un modèle dyna¬ 
mique qui décrit les fluctuations temporelles des variables du système. 


MODÈLE ÉCO-ÉNERGÉTIQUE 

Jusqu’à une époque relativement récente, les travaux traitant à la fois de l’écologie 
et de l’énergétique concernaient principalement l’étude des flux d’énergie le long des 
chaînes trophiques, pour des populations déterminées d’espèces animales dans des 
systèmes naturels. Les problèmes actuels d’approvisionnement en énergie de même 
que l’intérêt croissant porté aux systèmes agro-alimentaires ont amené le développement 
d’une nouvelle tendance qui étudie la totalité de la circulation de l’énergie dans les éco¬ 
systèmes modifiés par l’homme, notamment les agro-écosystèmes (p. ex. Pimentel et al. 
1973; Vieira da Silva & Souchon 1977). 

Une analyse éco-énergétique du système des étangs dombistes a été effectuée par 
Marchenay (1979). Ce modèle, dont la figure 2 présente l’application à l’étang 
Grand-Glareins, cherche à estimer le coût énergétique de la production alimentaire 
que fournit cet écosystème. L’analyse éco-énergétique s’efforce en effet d’identifier toutes 
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Fig. 2. 

Flux d’énergie dans l’agro-écosystème de l’étang Grand-Glareins. (1 an d’assec, 2 ans d’eau. 
Culture: avoine). Le bilan brut est calculé à partir des productions obtenues par l’exploitant. 
Le bilan net tient compte des transformations ultérieures nécessaires pour l’alimentation humaine. 
Ainsi, les céréales servant à l’alimentation du bétail sont transformées en lait et viande. 
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les sources et consommations d’énergie et leurs transferts le long des chaînes de pro¬ 
duction. 

Les informations concernant les entrées et sorties ont été obtenues par des enquêtes 
menées auprès des exploitants. Les données concernant le travail, la consommation 
de combustible, l’utilisation d’engrais et d’herbicides sont assez constantes d’année 
en année sur une période de dix ans, tandis que les productions sont plus irrégulières. 
Toutes les valeurs utilisées dans le modèle sont des moyennes sur 10 ans. 

Les entrées d’énergie dans le système sont multiples et pas toujours faciles à mesurer. 
Le soleil est à la base de tout écosystème. Il représente la source d’énergie la plus impor¬ 
tante, mais étant inépuisable et gratuit, il ne sera pas comptabilisé dans le coût énergé¬ 
tique des productions. 

Le travail humain est nécessaire à l’entretien et à l’exploitation du système (labours, 
moissons, pêches, entretien des fossés). L’énergie dépensée par l’homme pour ce travail 
est estimée par la fraction de la ration alimentaire effectivement transformée en travail. 
Elle est de 1254 kJ/h pour un travail agricole moyen (Passmore et al. 1974). 

Les machines agricoles nécessaires à l’exploitant fournissent de l’énergie méca¬ 
nique au système et consomment du combustible fossile pour leur fonctionnement 
et leur fabrication (Pimentel et al. 1973; Deléage et al. 1977). Les engrais et herbicides 
sont représentés dans le système par la dépense énergétique nécessaire à leur obtention 
(Marchenay 1979). 

Les sorties du système sont essentiellement constituées par des productions de 
poisson, céréales et paille. Leur contenu en énergie est estimé par leur valeur calorique, 
mesurée en bombe calorimétrique. 

Un budget « brut » peut être calculé à partir des sorties telles que l’exploitant les 
obtient sans les transformer. Le rendement brut sera le rapport entre la quantité 
d’énergie extraite et la quantité d’énergie introduite dans le système. L’énergie pro¬ 
duite et réintroduite sous forme d’alevins ou de semences appartient à la circulation 
interne du système et n’intervient pas dans le calcul du rendement. Les empoisonnages 
de carpes figurent comme une énergie introduite car elles sont élevées dans des étangs 
spéciaux, dits de pose ou d’alevinage, qui produisent les « panots » de deux étés. Dans 
l’exemple de l’étang Grand-Glareins, la production de poissons, après avoir retiré la 
fraction nécessaire au réempoissonnage s’élève à 1075 MJ/ha/an avec un rendement 
brut de 4,38. En assec, la culture de l’avoine, après prélèvement des semences, produit 
32’451 MJ/ha/an avec un rendement brut de 3,93. Le rendement sur un cycle complet, 
formé de deux ans d’évolage et d’un an d’assec est de 3,95 (sans la paille, laissée sur le 
terrain). 

Mais les étangs dombistes ayant d’abord une vocation alimentaire, il est plus 
intéressant de calculer un budget « net » sur la base de la production effectivement 
consommée par l’homme. En Dombes, les céréales subissent une nouvelle transformation 
comme nourriture pour le bétail, soit pour la production de lait, soit pour la production 
de viande. Par ailleurs, seule une partie du poisson est propre à la consommation humaine. 
Dans l’étang Grand-Glareins, la production nette de l’évolage est de 641 MJ/ha/an 
en poisson (partie comestible), avec un rendement net de 2,61. Lors de l’assec, la pro¬ 
duction est de 10T45 MJ/ha/an avec une option «lait», ou de 1368 MJ/ha/an avec 
une option « viande ». Les rendements nets de Passée sont de 0,56 pour le lait ou de 
0,11 pour la viande. Les rendements nets ainsi calculés restent approximatifs dans 
la mesure où le régime alimentaire du bétail ne peut être constitué uniquement des céréales 
et du fourrage ensilés. 

En ce qui concerne le rendement en protéines, les productions nettes sont de 
22 kg/ha/an sous forme de poisson comestible lors de l’évolage, et de 127 kg/ha/an 
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sous forme de lait, ou de 31 kg/ha/an sous forme de viande, lors de Passée. Les coûts 
énergétiques par kg de protéines produites sont respectivement de 11,4 MJ/kg pour 
le poisson, de 140 MJ/kg pour le lait et de 371 MJ/kg pour la viande. 

La production à l’hectare, en période d’assec, est inférieure de 10 à 20% à celle des 
terres avoisinantes. En contrepartie, le rendement énergétique de l’étang cultivé est 
plus intéressant. Les comparaisons manquent, mais on peut citer à titre indicatif quelques 
rendements énergétiques bruts, à comparer au chiffre de 3,93 obtenu pour l’avoine dans 
l’étang Grand-Glareins : 

Maïs-grain, Etats-Unis: 2,82 (Pimentel et al. 1973) 

Maïs-grain, Angleterre: 2,32 (Leach 1976) 

Maïs-grain, France: 3,20 (Mercier 1978) 

Avoine, France: 3,40 (Mercier 1978) 

L’évolage, qui nécessite un apport énergétique particulièrement faible (peu de tra¬ 
vail mécanique, pas d’engrais), présente un rendement énergétique très élevé. Nous 
avons vu que l’étang Grand-Glareins assure un rendement brut de 4,38 et un rendement 
net de 2,61. Il est spectaculaire de comparer ces chiffres à ceux obtenus par Leach (1976) 
ou Angelidis (1978) pour la pêche en mer: respectivement 0,05 et 0,07. 

L’analyse éco-énergétique permet de comparer différents agro-écosystèmes du 
point de vue de leur rendement énergétique. Elle fournit ainsi une base de comparaison 
indépendante des problèmes économiques et monétaires spécifiques à chaque pays, 
ce qui la rend utilisable à très large échelle. Mais elle ne peut bien s’appliquer qu’à des 
systèmes à l’état stationnaire et n’a qu’une faible valeur prédictive. Elle ne peut répondre 
à des questions concernant la dynamique du système comme, par exemple, ce qui se 
passerait si le cycle évolage-assec était modifié. Le second modèle présenté dans ce 
travail tente de fournir un outil prédictif. 


MODÈLE DYNAMIQUE 

Comme nous disposions de fort peu de données quantitatives relatives aux variables 
biologiques des étangs, nous avons décidé de n’obtenir aucun paramètre par calibration. 
Ceci nous a contraints à élaborer un modèle suffisamment détaillé pour pouvoir y 
incorporer des constituants connus par des études sur d’autres systèmes. Il existe en 
effet une abondante littérature sur la modélisation des écosystèmes aquatiques, parmi 
laquelle on peut citer les travaux bien documentés de Chen & Orlob (1975), Di Toro 
et al. (1977), Jdrgensen (1976, 1980), Jdrgensen et al. (1978), Nyholm (1978), Scavia 
(1980), Scavia & Park (1976), Thomann et al. (1975). Une comparaison de nombreux 
modèles a été faite par Swartzman (1979) et Swartzman & Bentley (1979). 

Le modèle comprend six compartiments (Fig. 3): phytoplancton, zooplancton, 
poisson, détritus, sédiment et azote minéral. Nous n’avons considéré que ce nutriment 
du phytoplancton car c’est le seul qui soit limitant dans le cas des étangs de la Dombes 
(Balvay 1980). Les six variables sont exprimées en milligrammes d’azote par litre d’eau. 
Pour des raisons d’homogénéité, cette unité est également adoptée pour le sédiment. 
Sa concentration réelle au fond, par unité de surface, s’obtient facilement en multi¬ 
pliant par la profondeur de l’étang. 

Il est fait l’hypothèse que l’ensemble des espèces constituant le phytoplancton sont 
consommées par le zooplancton. En effet, l’apparition d’espèces non consommables, 
telles les fleurs d’eau, est un phénomène rare que les exploitants s’efforcent d éviter. 

I 
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De même, l’ensemble du zooplancton est supposé être consommable par les poissons. 
Inversement, ces derniers sont supposés être uniquement planctivores, les poissons 
prédateurs (brochets) étant négligés. Le « détritus » comprend le détritus organique 
en suspension et la couche superficielle du sédiment. 



Fig. 3. 

Modèle du cycle de l’azote. Les six compartiments sont: phytoplancton, zooplancton, 
poisson, détritus, sédiment et azote minéral. Dans les étangs de la Dombes, 
l’azote minéral est le seul nutriment limitant du phytoplancton. 


La dynamique du modèle est régie par le système de six équations différentielles 
du premier ordre illustrées par la figure 3. La signification des flux représentés par 
les flèches est la suivante. PRODP est la production primaire nette. Il n’y a donc pas de 
respiration du phytoplancton. ZORES et PORES sont les respirations du zooplancton 
et du poisson. PHMOR et ZOMOR sont les mortalités. PHMOR comprend également 
la sédimentation du phytoplancton. Aucune mortalité n’est prévue pour le poisson car, 
normalement, celui-ci ne fait que croître au cours de l’année d’exploitation, à l’issue 
de laquelle il est intégralement pêché. ZOPRD et POPRD sont les prédations, ZOASS 
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PRODUCTION PRIMAIRE 

PRODP = PPM4X*pplum*ppazo*pptmp*phy 

1 * RAY(t) + RAYK 


Où PPLUM = 


PROF* EXT 


£n 


RAY(t) e PP0P * EXT + raYK 


EXT = EXT/) + EXT1*PHY 
AZO 


PPAZO 


AZO + AZOK 


MORTALITES 

PHMOR = PHMORT* PMTMP*PHY 
ZOMOR = ZOMORT *ZMTMP * ZOO 

RESPIRATIONS 

ZORES = Z ORESP* Z RTMP * ZOO 
PORES = P0PPSP*PRTMP*POI 


PREDATIONS 

PH Y 

ZOPRD = ZPMAX *-—— * ZPTMP * ZOO 

PHY + PHYK 

ZOO 

POPRD = PZMAX * — --* PZTMP * POI 

ZOO Z DU K 

SEDIMENTATION 

DTS' = ( DETSED/PROF ) * DET 

MINERALISATIONS 

DTMIN = DET MI R *DMTMP*DET 

SDMIN = SEDMIN *SMTMP*DTSED 


INFLUENCE DE LA TEMPERATURE (loi de Van't Hoff) 
(TMP(t)- 20) 


xxTMP = PPxa: 


(xx = PP,PM,ZM,ZR,PR,ZP,PZ,DM,SM) 


Tableau 1. 

Expression détaillée du modèle. Les quantités en italiques sont les paramètres et les variables 

exogènes. 
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et POASS étant les taux d’assimilation. Enfin, les flux DTSED, DTMIN et SDMIN 
représentent la sédimentation du détritus, la minéralisation du détritus et celle du sédi¬ 
ment. 

Le tableau 1 donne l’expression détaillée de tous les flux, les quantités en italiques 
étant les paramètres et les variables exogènes. Nous avons adopté des interactions de type 
multiplicatif lorsque plusieurs facteurs agissent en même temps. La température agit 
selon la loi de Van’t Hoff sur tous les flux, à l’exception de la sédimentation du détritus. 

Les consommations (prédations du poisson et du zooplancton, absorption d’azote 
par le phytoplancton) sont supposées répondre à la cinétique de Michaelis-Menten qui 
limite la quantité consommée f(x) à une valeur de saturation 

(1) 


où x est la densité de l’aliment et xk la densité de demi-saturation. Cette loi est celle 
qui est adoptée par la quasi-totalité des modèles de simulation des systèmes hydrobio¬ 
logiques. 


Cette même loi est adoptée pour l’influence du rayonnement lumineux sur la 
photosynthèse (loi de Baly) 


/(/) = 


/ 

/ + RAYK‘ 


( 2 ) 


Mais l’intensité lumineuse I s’atténue avec la profondeur selon la loi de Beer 

I(z) = RAYe- EXT - z , (3) 


où RAY est l’intensité du rayonnement transmis à la surface de l’eau, et EXT est le 
coefficient d’extinction. Ce dernier est supposé recevoir deux contributions; l’une cons¬ 
tante, due à l’eau, et l’autre proportionnelle à la concentration de phytoplancton (auto¬ 
ombrage). 

En supposant que le phytoplancton soit distribué verticalement de façon uniforme, 
il est nécessaire de faire une moyenne de (2) en calculant l’expression 


PPLUM = 


1 

PROF 


dont la forme explicite figure au tableau 1. 


PROF 
fV(z))dz , 
0 


(4) 


La loi de Baly est assez largement utilisée, mais certains lui préfèrent la loi de 
Steele (1962) qui décrit une chute de la productivité au-delà d’une intensité lumineuse 
optimale. Bien que ce phénomène soit exact pour chaque espèce d’algue prise séparé¬ 
ment, il est vraisemblable que la composition du phytoplancton soit telle qu’il se trouve 
toujours une espèce pour exploiter au mieux le rayonnement disponible, la réponse 
résultante ne montrant alors plus de chute de la productivité. 

La minéralisation du sédiment est d’une structure différente des autres processus. 
L’expression adoptée a été proposée par Nyholm (1978) pour rendre compte du fait 
que ce taux de minéralisation est positivement corrélé à la température et à la sédimen¬ 
tation. Ceci permet d’éviter de devoir décrire de façon détaillée les processus internes 
au sédiment. 
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Paramètre 

Valeur 

Source 

PROF 

1,2 m 

„ r . -1 

estimation moyenne 

PPMAX 

2,5 jour 

-2 -1 

Thomann et al. 1975, Di Toro et al. 1977, 

Scavia 1980, Handbook A57,A131 

RAYK 

1003 J cm jour 

-1 

Tiwari 1978 

EXT0 

0,27 m 

-1 -1 

Jurgensen 1980 

EXT1 

6,0 rag m 

Jurgensen 1980 

PHMORT 

0,1 jour 1 

Handbook A72 

ZOMORT 

0,005 jour 

-1 

Handbook A458 

ZORESP 

0,02 jour 

Handbook A471, Chen & Orlob 1975 


„ -1 


PORESP 

0,002 jour 

-1 

Handbook A222 

DETMIN 

0,2 jour 

Handbook G104, Jurgensen 1980 

SEDMIN 

1/0 

Nyholm 1978 

DETSED 

0,003 m jour 
~ -1 

Jurgensen et al. 1978, Handbook A202 

ZPMAX 

0,36 jour 

„ , -1 

Jurgensen et al. 1978, Handbook A440 

P Z MAX 

0,14 jour 

-1 

d'après Handbook A215 

AZOK 

0,2 mg 1 
„ n -1 

Jurgensen et al. 1978, Handbook A48 

PHYK 

0,2 mg 1 

-1 

Jurgensen et al. 1978, Swartzman & Bentley 1979 

ZOOK 

0,35 mg 1 

J0rgensè.: et al. 1978 

ZOASS 

0,6 

Handbook A203 

POASS 

0,4 

Handbook A203 

THPP 

1,082 

Handbook Al31 

THSM 

THPM 

THZM 

THZR 

1,14 

Nyholm 1978 

THPR 

THZP 

THPZ 

THDM 

1,045 

Di Toro et al. 1977 

Teneurs en 

azote 


phytopl. 

7,0 N/chloro a 

Handbook A75 

phytopl. 

0,08 N/poids sec 

Handbook A74 

zoopl. 

0,008 N/poids sec 

Handbook A459 

poisson 

0,11 N/poids sec 

Handbook A218 

poisson 

0,026 N/pds humide 

Handbook A217 


Tableau 2. 

Valeurs des paramètres. 
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SIMULATIONS 

Pour la plupart des paramètres, nous ne disposions pas de valeurs mesurées sur le 
terrain. C’est pourquoi il a fallu adopter des estimations tirées de la littérature, en compa¬ 
rant de nombreuses sources. Certains paramètres peuvent varier de façon relativement 
importante selon les auteurs et nous avons dû quelquefois faire un choix, nous effor¬ 
çant d’adopter les valeurs les plus largement utilisées. Les valeurs finalement adoptées 
pour les simulations sont énumérées au tableau 2, où seules figurent les sources les 
plus fréquentes. Le « Handbook » dont il est souvent fait mention est un gros ouvrage 


Année 

1972 

1973 

1976 

1977 

1978 

1979 

1980 

Empoissonnage (kg) 

1249 

1070 

930 

955 

898 

873 

948 

Empoissonnage (mgN/1) 

0,11 

0,09 

0,07 

CD 

O 

O 

0,07 

0,07 

o 

o 

00 

Pêche (kg) 

4074 

2749 

5904 

5494 

6550 

4559 

2870 

Pêche (mgN/1) 

0,32 

0,22 

0,47 

0,44 

0,52 

0,36 

0,23 

Lisier (t) 

29,7 

29,7 

29,7 

29,7 

29,7 

29,7 

29,7 

Lisier (mgN/1) 

0,16 

0,16 

0,16 

0,16 

0,16. 

0,16 

0,16 


Tableau 3. 

Données relatives à l’étang du Grand Guerrier (27 ha de superficie, 1,2 m de profondeur moyenne). 


édité par Jurgensen (1979), qui rassemble un très grand nombre de paramètres recueillis 
dans la littérature. Comme toutes les concentrations ont été exprimées en azote, il a 
été nécessaire d’utiliser les facteurs de transformation mentionnés au tableau 2 lorsque 
les sources utilisaient d’autres unités. 

A titre d’exemple, nous avons choisi de simuler l’étang du Grand Guerrier sur la 
période 1972-1980, les années 1974 et 1975 étant des années d’assec. Le tableau 3 présente 
les données connues sur cet étang. Les paramètres ayant été fixés (tableau 2), il restait 
à connaître les variables exogènes et les valeurs initiales des six variables d’état. Concer¬ 
nant les premières, nous disposions des données suivantes: quantité de lisier déversé 
dans l’étang, températures moyennes à Marlieux et rayonnement solaire global à Mâcon. 

Les valeurs initiales du poisson sont évidemment celles de l’empoissonnage. Pour 
les autres variables, il se posait le problème dû au fait qu’un étang vidé lors de la pêche 
est rempli par l’eau provenant d’un autre étang ou du bassin versant. Pour simplifier, 
nous avons pris le cas fictif où l’étang serait rempli par sa propre eau (ou, ce qui revient 
au même, serait pêché sans être vidé). A l’exception de la première année, les valeurs 
initiales seront donc déterminées par les valeurs finales de l’année précédente. 

La valeur initiale du détritus doit être augmentée de la quantité de lisier de vache 
déversé dans l’étang (1,1 t/ha), dont la teneur en azote est en moyenne de 1,7 g/kg 
(G. Balvay, comm. pers.). Cela représente un apport de 0,16 mgN/1. 
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Pour 1972, première année de la simulation, nous avons pris des valeurs initiales 
jugées moyennes: 0,3 pour le sédiment, 0,2 pour l’azote, 0,3 pour le détritus, 0,1 pour le 
phytoplancton et 0,0005 pour le zooplancton. En 1976, à la remise en eau après Passée, 
nous avons fait l’hypothèse que l’étang était rempli par une eau identique à celle de la 
fin de 1973, à la différence que tout le sédiment est supposé s’être transformé en azote 
minéral. 

Un scénario de valeurs initiales ayant ainsi été fixé pour toutes les années, les simu¬ 
lations sont effectuées séparément pour chaque année avec un pas journalier. La 
figure 4 donne l’exemple de la première année. 


1972 



JOURS 


Fig. 4. 

Simulation de l’étang du Grand-Guerrier pour l’année 1972. 

La biomasse de poisson au 365 e jour est, par définition, la « pêche prédite ». 


DISCUSSION 

La figure 5 présente la comparaison des valeurs annuelles finales du poisson, et les 
pêches effectives dans l’étang du Grand-Guerrier (sans le brochet). Compte tenu du 
nombre d’hypothèses simplificatrices qu’il a fallu faire, il ne fallait pas s’attendre à un 
accord parfait. Dans ces conditions, il est déjà satisfaisant que les pêches prédites 
soient du même ordre de grandeur que les pêches réelles. De plus, le modèle rend compte 
d’un important phénomène: la baisse graduelle des rendements tout au long de l’évolage, 
et un bond très net après l’assec. Ceci est dû au blocage des nutriments dans le sédiment, 
et à leur libération lors de la mise en culture. 

Il resterait beaucoup à faire pour aboutir à un modèle robuste et fiable. Rappelons 
qu’aucun paramètre n’a été fixé par calibration. La plupart d'entre eux ont été tirés 
de la littérature et ne correspondent peut-être pas aux valeurs du système étudié. Des 
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estimations directes du plus grand nombre possible de paramètres seront les bienvenues. 
En particulier, ceux relatifs à la prédation du poisson sont mal connus. Il serait alors 
possible de décrire plus finement le compartiment « poisson », en introduisant par exemple 
le brochet. 

Il conviendrait aussi de donner de Tassée une description moins grossière que celle 
qui a été adoptée, à savoir une simple minéralisation du sédiment. L’assec s’accompagne 



ANNEE 

Fig. 5. 

Etang du Grand-Guerrier. L’accord entre pêches réelles et pêches prédites n’est pas parfait, 
mais le modèle rend compte d’un phénomène important: la baisse graduelle des rendements 
tout au long de Tévolage, et un très net bond après î’assec. Ceci est dû au blocage des nutriments 
dans le sédiment, et à leur libération lors de la mise en culture. 


également d’un apport d’engrais et d’une production de détritus. Si l’on obtient en 
outre une estimation de la production de l’assec, il serait alors possible d’utiliser le 
modèle pour trouver la politique d’assolement qui maximiserait la production totale. 

Ainsi, par son caractère prédictif, l’approche en simulation dynamique joue un rôle 
différent mais complémentaire de celui de l’approche éco-énergétique qui, elle, donne 
une description essentiellement statique. Par là-même, l’analyse éco-énergétique est 
relativement facile à effectuer puisqu’il suffit de mesurer les flux et niveaux d’énergie 
et de matière, tels qu’ils existent. L’approche dynamique nécessite en plus la connais¬ 
sance des lois auxquelles répondent les divers composants du système lorsqu’ils sont 
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placés dans des situations différentes de celles qui existent. Cette difficulté conduit 
presque toujours à élaborer un modèle plus simplifié que lors de l’approche éco-éner- 
gétique. 

Cette dernière permet en revanche de tenir compte d’un grand nombre de compo¬ 
sants et donne une vision très complète de la circulation de l’énergie. Dans le cas de 
la Dombes, elle a montré la grande efficacité du système. Associée à une exploitation 
agricole, la culture des étangs tend à faire augmenter les rendements énergétiques glo¬ 
baux de manière très sensible, en transformant rapidement une productivité primaire 
et secondaire élevée en produits alimentaires de valeur pour l’homme. Dans un milieu 
où la densité humaine est relativement faible, on a là un type d’exploitation extensive 
très intéressant. 

Deux problèmes d’importance vitale vont se poser de façon aiguë au cours des 
prochaines années : 

— Le caractère renouvelable des entrées d’énergie est certainement l’un des plus 
importants de l’agriculture. L’exemple des étangs nous montre qu’il existe des solu¬ 
tions qui, sans être un retour à une technique révolue, offrent une possibilité de gérer 
de façon optimale une entreprise basée sur l’exploitation d’une productivité naturelle. 

— La désertification des sols est liée aux problèmes énergétiques. Si l’on doit dimi¬ 
nuer les apports d’engrais et les coûts énergétiques de combustibles fossiles, il faut trou¬ 
ver des solutions d’urgence pour maintenir la fertilité des terroirs. L’analyse éco-éner- 
gétique peut y contribuer efficacement, en même temps que des études poussées sur l’im¬ 
portance de l’humus, sur les cycles bio-géochimiques et notamment celui de l’azote. 


Résumé 

La Dombes est un plateau situé au sud du département de l’Ain, dont le centre 
est une région d’étangs artificiels, exploités alternativement en eau et en assec. L’exploi¬ 
tation en eau (évolage), qui dure deux à six ans, est consacrée à la pisciculture de carpes, 
blancs et brochets. En principe, les étangs ne sont pas nourris artificiellement. La 
pêche annuelle se fait par vidange de l’étang. L’exploitation en assec, qui dure un ou 
deux ans, est consacrée à la culture de l’avoine, du blé ou du maïs. 

Deux types de modèles de ce système sont présentés dans ce travail. Le premier est 
un modèle éco-énergétique qui décrit le fonctionnement du système dans un état station¬ 
naire, sous l’angle particulier des flux d’énergie. Il permet d’estimer le coût énergétique 
de la production alimentaire que fournit cet écosystème. 

Le second modèle est de type dynamique. Il décrit les fluctuations temporelles des 
variables du système, lors de l’évolage. Il est constitué des six compartiments phytoplanc- 
ton, zooplancton, poisson, détritus, sédiment et azote minéral. La dynamique de ces 
six variables d’état est régie par un système d’équations différentielles. La forme ana¬ 
lytique des interactions, ainsi que les valeurs des paramètres, reposent sur des données 
de la littérature. Après simulation de sept années, il est possible de comparer les pro¬ 
ductions prédites par le modèle et les productions réelles. L’accord n est pas parfait 
mais il est assez bon qualitativement. Le modèle rend compte de la baisse graduelle de 
rendement au cours de l’évolage. 
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